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^ ■ Abstract 

(N 

> , 

f\J We discuss the implications of CP violation as well as final state interaction 

\ phases in the experimental search for — mixing. At the present level of 

o 



sensitivity, these are not yet a significant systematic experimental limitation. 



Oh! I. INTRODUCTION 



- ^ , As was recently noted by Blaylock, Seiden, and Nir [|l| due to final state interaction (FSI) 

^ , a term proportional to AM t e may appear in the rate of wrong sign D decays even in 
the absence of CP violation. Moreover, in some extensions of Standard Model which have 
large values of both AM and significant CP violation, a similar term may arise. Blaylock 
et al. have suggested that a value of AM larger than the present experimental limit can 
be accomodated if one of these previously neglected terms destructively interferes with the 
other time dependent terms which arise from mixing (proportional to e~'"*) and from 
doubly Cabibbo suppressed decays (DCSD) (proportional to e"^*). They suggest that this 
may invalidate the use of existing limits from time dependent mixing studies at fixed target 
experiments [0 , P] to constrain extensions of the Standard Model. 
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Below, we give expressions for the time dependence in the general case and then attempt 
to estimate the maximum size of the terms proportional to te^^*. 



II. FORMALISM FOR MIXING 

We follow the notation of references ||^, 0. Let the mass eigenstates be Ds, Dl- 
Then 

\Ds >=p\D° > +q\D^ > 



\Dl >= p\D'^ > -q\D° > 

In the limit of no CP violation, p=q=l/ 

Let AM = Ml — Ms and AF = F^ — r5 denote the mass difference and lifetime difference, 

respectively. Let A denote the amplitude for < f\H\D'^ >, B the amplitude for < f\H\D^ >. 

Let A = -— and A = The decay rate is then given by 
q B p B 

p-^* o AF^ 
r{D\t) K+7T-) = ^\B\^\^^{A\X\^ + (AM^ + -^)t^ + 2Re{X)Art + 4Im{X)AMt} 

(1) 

up to terms of order The decay rate for the charge conjugate reaction is given by the 

same expression replacing A with A, B with B, and q/p by p/q. 

p"^* - n - AF^ 
r{D\t) K-7f+) = l_|5n^:|2{4|A|2 + (AM^ + =^)f + 2i?e(A)AFt + 4Im{X)AMt} 



(2) 



In the past, it was assumed that the term proportional to AM t changes sign when averaging 
over a sample with equal numbers of and mesons 0, 0. This assumption is not 
correct in general as was noted in Reference 

The previous experimental analyses 0, 101) lH considered the deterioration of the limit 
in the case when the term proportional to AF t interfered destructively with the mixing and 



DCSD components. The Standard Model expectation for AF is many orders of magnitude 
below the current experimental sensitivity so this interference scenario is very unlikely. In 
most new physics scenarios which would give Vmix ~ 0(10^^), AF is not enhanced whereas 
values of AM much larger than those expected from the Standard Model are possible. It is 
also possible to experimentally verify that AF can be neglected by measuring the D meson 
lifetime in a CP eigenstate e.g. —>■ K~K^ and comparing to the lifetime in K^tt~^ 

i, i- 

We now consider equations (||), (H) in the following situation. Let 

Q 

- = ae^' 
B 

and = — ; — -; r. The phase ch is due to CP violation in the mass matrix. A non- 

F(DO ^ K-7r+) ^ ^ 

zero value of 6 may arise if the amplitudes A and B have different FSI. Alternatively, if there 

are complex contributions to one of the amplitudes (e.g. A) that are not present in the other 

(e.g. B), this can lead to an overall phase in — . We have assumed that there is no direct 

B 

CP violation in the amplitudes (and hence e.g. r{D^ K^tt^) = r{D^ K^tt^)) In 
addition, we neglect the small phase in A/B from the CKM matrix, which is approximately 
A'^X'^rj in the Wolfenstein parameterization and which lies in the range (2.3 — 5.3) x 10~^. 
The decay rate for wrong sign decays to K^tt' is given by 

F(D°(t) ^ K+TT-) = X 

^^{4a'p' + {AM' + ^)t' 

+2ap cos(0 + 5) {Art) + 4 sin(0 + 5)apAMt} 

The corresponding rate for the charge conjugate reaction is obtained by replacing the 
phase from CP violation with —0 and by changing (3 to 



r{D''{t) K-71+) = ^\B\' X 



+2^ cos(-(/) + 5) (AH) + 4sm(-0 + 5)^/^Mt} 

In the experimental analyses, the time dependent rate integrated over both types of 
particles is used: 

r(D°(t) ^ K+TT-) + T{D%t) ^ K-TT^) 
This rate, which will be denoted by r(D^(t) + D^(t)), is given by 

r{D\t) + D\t))^2F{t)x (3) 

+q;(/3cos(-0 + (^) + 4 cos{(f) + 5))Ar t 

P 

+2a{psm{-(p + S) + ^sm{(p + 6))AM t} 

where F{t) = ^e"^* 



III. EFFECTS OF FSI AND CP VIOLATION 

Two scenarios are considered in what follows. First the case of no CP violation but 
significant final state interactions (FSI) and then the case of both large CP violation and 
significant final state interaction are examined. 



A. Effects of FSI 

In the first scenario, consider the case of large mixing with AM >> AMsm, the value 
in the Standard Model. Assume that this does not lead to an enhancement of AF i.e. 
AF^M = AF << AM and allow for non-zero 5 but no CP violation (0 = 0,/3 = 1). The 
above equation then reduces to 

AF^ 

V{D\t) + D\t)) = 2F{t){W + (AM^ + + 4Q;(sin(5))AM (4) 
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In order to determine the size of the new term proportional to AM t, the values of 
the phase difference 6 are considered in various models. This will allow an estimate of 
the additional experimental systematic error that is incurred from ignoring FSI. This phase 
difference 6 is zero in the limit of exact SU(3) symmetry. The values of 6 from various models 
are given in Table |. Large values of the phase 6 occur when SU(3) breaking is largest. We 
use the experimental result from CLEO II for K+n- (a^ = 0.0077 ± 0.0025 ± 0.0025) 

and assume that it is entirely due to DCSD. This is found numerically to give the most 
conservative upper limit on the size of the interference effect. 

In general, the amplitudes for the K^tc^ and K^tt^ can be written as: 



where Ai and A^ are the quark decay contributions into J = 1/2 and 7 = 3/2 final states 
respectively. C is the W-exchange contribution and 6i and ^3 are the FSI phases. Ai, C are 
the corresponding DCSD amplitudes after the CKM factor —6l has been factored out. The 
phase shifts in a given isospin eigenstate for particles and antiparticles are identical by CPT 
invariance (which we assume as stated explicitly). 

The phase 5 vanishes if two conditions are satisfied: (i) (5i — ^3 = 5i — ^3 and (ii) 
A^/Ai = A3/A1. The first condition follows from CPT invariance and the second is satisfied 
if SU (3) symmetry holds. Hence if SU (3) is an approximate symmetry, the phase 6 should 
be small. The models used have been tuned to reproduce the observed magnitude of SU (3) 
breaking in D decays. To obtain more information, we turn to the detailed model fits. 



A{D^ K'7T+) = e'^'[{Ai + C)e*('^^-^^) + A3] 



(5) 



^ K+71-) = -ey^'[{Ai + (7)6^(^1-^3) + ^3] 



(6) 



B. Details of the Models 



In the model of Chau and Cheng, 



and 



5i-53^ 90°, ^3 ^ 0. (9) 

Then 

A{D'> ^ K-7r+) ^ {0.72)e' (10) 

A{D^ ^ K+TT-) ^ -el{im)e' (11) 

This yields a phase difference between the two decay modes of 5 = 4°. If the W-exchange 
contribution C is omitted, the phase difference becomes 6 — 5°. 
In the model of Buccella et al., one has 

Ai ^ 4.35, A^ ^ -2.3, C ^ -0.5 (12) 

ii ^ 5.2, As = -2.3, C^-C (13) 

and 

5i-53^ 25°, ^3 0. (14) 

Then 

A(L>° ^ X-7r+) ^ (2)e' ^^-'^ (15) 

>1(L'° ^ i^+TT") ^ -^,'(3.7)e^ (16) 

leading to a phase difference of 5 = 13°. Omitting the W-exchange term gives a slightly 
smaller value of 6°. It should be noted that a, 5i — 25° relates to Sr for the 7 = 1/2 0+ 
resonance in the Kn channel by 
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where B = BR{0+ Kn) ^ 0.50, T ^ 200 MeV, A = Mr - ^ 70 MeV and 6r = 
(55 - 65)°. 

The models discussed above predict 

^ ^ - (2.3-3.4)tan^^, 



c 



which is compatible with the CLEO II measurement. There are also other models for D 
decays in which a value for the phase difference S can be extracted WM. Since it is difficult 



to assign errors to these predictions, we regard 0° — 13° as a reasonable range for 6. In 
order to explore the range of 6 in the models, we have calculated the value of S omitting the 
W-exchange term. This corresponds to a dramatic change in the parameters of the models. 
TABLE I. Values of 6 in various phenomenological models of D meson decay. 



Exact SU(3) limit |T| 0° 

Chau and Cheng |l| 4° 

Chau and Cheng (no W-exchange) ||l^ 5^^ 

Buccella et al. 13° 

Buccella et al. (no W-exchange) 6° 



C. Summary of the Interference Effect from FSI 

To summarize, the phenomenological models which have been tuned to agree with the 
observed branching fractions and the fits to the D meson data give 6 in the range of 5° — 13°. 
To evaluate the possible experimental consequences, consider the case of maximal destructive 
interference (0 = 0, /5 = 1), with 6 = 13°. We allow a one standard deviation variation on 
Rdcsd = from the CLEO II measurement in order to obtain an upper limit on the effect 
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of the interference term. We set r^ix, the ratio of integrated rates for mixed events relative 
to unmixed events, to the E691 upper bound [|T^. The contributions of the mixing term, 
the DCSD term, and the term proportional to AMt are shown in Figure 1. These time 
dependent searches are most sensitive to excess events from mixing for t > 0.22 ps = 
where the combinatorial backgrounds are manageable and where the mixing term is expected 
to peak. In addition, there is no loss in efficiency for the mixing component when this cut 
is imposed. An upper limit of t < 4.0 ps is also imposed. The change in the observed event 
yield for various values of Rdcsd and maximal destructive interference are given in Table U 
These were calculated for the scenario with maximal destructive interference and 6 = 13°. 
We also give the change in the observed event yield for t > 2t£)0 (this is the region where 
mixing peaks and the experiments are most sensitive) in Table |T|. This change is at most 
10-15% and is well within the experimental systematic error assigned by the E691 and E791 
experiments to their limits. 

TABLE II. The change in wrong sign event yield for t > 0.22 ps with maximal destructive 

interference, Vmix = 0.37%, and 6 = 13". 

Rdcsd a Yield (%) 

0.0052 10% 

0.0077 8% 

0.0102 1% 

TABLE III. The change in wrong sign event yield for t > 0.88 ps (2t£,o) with maximal destruc- 
tive interference, rmix = 0.37%, and 5 = 13°. 

Rdcsd A Yield 

0.0052 12% 
0.0077 10% 
0.0102 9% 
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IV. EFFECTS OF CP VIOLATION 



1 



Now consider the contribution of CP violation. Let (3 = 1 — e and — = 1 + e [|T6[. We 

P 

assume e is small compared to 1 and retain only terms linear in e; this is justified in the SM 
and even more so when AM is enhanced and AF/AM << 1. We allow the phase to be 
arbitrary. With these definitions and AF << AM, the expression for T[D^[t) + D^{t)) now 
becomes: 

F(D0(t) + 50(t)) = (19) 

AF^ 

2F{t) X + (AM2 + 

+a(cos(-0 + S)+ cos(0 + 5))AF t + ae(cos(0 + 5)- cos(-0 + 5))AF t 
+2a(sin(-0 + 5) + sin(0 + 5))AM t + 2ae(sin(-0 + 5) - sin(0 + (5))AM t} 

The quantity e is assumed to be small as in Ref however, the CP violating phase 
can be large as is the case for certain extensions of the Standard Model. The quantity e for 
D mixing is given by |[T^ 

^ 2 ^ (20) 



1 AM2 + 



in the Standard Model and is already small (e < 0(2%)) |T^. 
In new physics scenarios with AF^m = << AM, 

(21) 

^ AM ^ ' 

For non standard models with Im(Mi2)/AM of order unity, tan0 may be large (0(1)). By 
contrast, 

, AF , Im(Mi2) , AF , 

e^2 — L^~2 «1 22 

^AM^ AM ^AM^ ^ ' 

The crucial point is that e is proportional to 1/AM and is highly suppressed if tan0 is of 
order unity and AM is enhanced. It is important to note that while tan(0) can be much 
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larger than the Standard Model expectation e will be even smaller than the value in the 
Standard Model for new physics scenarios in which AM is enhanced. 
The total wrong sign rate can then be reduced to 

r{D\t) + D\t)) = 2F{t) (23) 
{W + (AM^ + ^)t'' 

+2a(cos(0) cos5))Ar t - 2ae(sin(0) sin(5))Ar t 
+Aae{cos(5) sin(0)) AM t + 4Q;(sin(5) cos(0)) AM t} 

With e as given above and AF << AM, the expression for r{D^{t) + D^{t)) becomes 

r{D\t) + D\t)) = 2F{t) 

+Aa{sm{S) cos(0))AM t} 

Hence, the term due to CP violation is too small to be observable when AM and Im(Mi2) 
are enhanced. 

As experimental sensitivity improves and become sensitive to mixing at the level Vmix < 
10"^, it is possible that better sensitivity to — mixing can be achieved by fitting the 
time distribution of r{D^ K'^tt~) — r(_D° — > K~tt'^). This rate, which will henceforth be 
denoted T{D^{t) - -D°(t)), is given by 

r(L>°(i) - L>°(t)) = 2F{t) X (24) 
{2Q;e(cos(0) cos5))Ar t - 2a{sm{(f)) sin(5))Ar t 
+4a(cos(5) sin(0))AM t + 4Q;e(sin(5) cos(0))AM t} 

In the limit that AF << AM and is large, this reduces to 

T{D%t) - D%t)) ^ 2F(t){4a(cos(5) sin(0))AM t + 4ae(sin(5) sin(0))AM t} (25) 

or neglecting the small term proportional to esin(5), 

10 



T{D%t) - D%t)) ^ 2F{t) [4a(cos((5) sin(0)) AM t] (26) 

Note that in this case, the long hved tail of DCSD does not contribute to the signal. In 
addition, as noted by Wolfenstein |T^, for small values of AM, the term proportional to 
AM t will be larger than the term in T{D^{t) + D^(t)) which is proportional to (AM t)^. 
This feature is illustrated in Figs. 2 (a), 2 (b). 



V. CONCLUSIONS 

The formalism presented here must be modified for the case of multibody modes such as 
D° fC+TT^TT^ or K^7i^7i^7T~ . FoT thcsc other modes, an additional complication 

is that the value of the final state phase difference, 6, may be different from the value in 
the case of K^n^/D^ K^-k^ and is not guaranteed to be small. It should also be 

remembered that limits on — mixing from studies of semileptonic decays do not have 
the complications from DCSD and other hadronic effects discussed here. 

At the present level of sensitivity and with reasonable (though model dependent) values 
for the phase difference 5, the AM t term which arises from FSI does not dramatically change 
the observed event yield for experiments which study the time dependence of mixing and is 
not yet a significant systematic experimental limitation. We suggest that future experiments 
determine systematic errors on their limits by using an upper limit on the phase difference 
5. 

The contribution from the corresponding term proportional to AM t due to CP viola- 
tion which arises in extensions of Standard Model is highly suppressed. This term is not 
observable at the present level of experimental sensitivity. However, as emphasized by Liu 
1^ and by Wolfenstein |T^, this term should not be neglected as experimental examination 



of the D^{t) - D^{t) distribution may allow more sensitive searches for — mixing in 
the future if the CP violating phase is large. 
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FIG. 1. The case of maximal destructive interference. The time dependence of DCSD(open 
circles), mixing(dotted), interference(dash-dotted) mixing and DCSD without interference (solid 
points), mixing and DCSD with interference (dashed) are shown. For this plot, the mixing rate is 
taken to be AM/F = 0.086 {vmix = 0.37%) and the DCSD rate is taken to be the central value 
of the branching fraction for D° K+Tr~ determined by CLEO II (i.e. Rdcsd = 0.0077). Other 
scenarios are discussed in the text. 
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FIG. 2. The time dependence for wrong sign events for the case of small mixing (AM/F = 0.02) 
for (a) T{D^{t) + where the dotted component is the FSI interference term and the dashed 

component is the usual mixing term (both scaled up by a factor of two in order to be visible), and 
for (b) F(D°(t) — 5°(t)). Note that in (b) there is no background to mixing from DCSD. 
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